The New Eclipsing Cataclysmic Variable SDSS 154453+2553 



m 
in 



Julie N. Skinner and John R. Thorstensen 

Dept. of Physics and Astronomy, 6127 Wilder Laboratory, Dartmouth College, Hanover, NH 

03755-3528; 
jns@dartmouth. edu 

Eve Armstrong 

Dept. of Astronomy, Columbia University, 550 West 120th Street, New York , NY, 10027 

and 
Steve Brady 

AAVSO, 25 Birch Street, Cambridge, MA, 02138 
ABSTRACT 



o 

(N 
l> 

o 

^ ■ The cataclysmic variable SDSS 154453+2553 was recently identified in the Sloan 

Digital Sky Survey. We obtained spectra and photometry at the MDM Observatory, 
which revealed an eclipse with a 6.03 hour period. The Ha emission line exhibits a 
strong rotational disturbance during eclipse, indicating that it arises in an accretion 
disk. A contribution from an M-type companion is also observed. Time-series photom- 
etry during eclipse gives an ephemeris of 2454878.0062(15) + 0.251282(2)E. We present 
spectroscopy through the orbit and eclipse photometry. Our analysis of the secondary 
star indicates a distance of 800 ± 180 pc. 



•rH . 1. Introduction 

X 

■ Cataclysmic variable stars (CVs) are close binary systems containing a white dwarf accreting 

matter from an extended secondary star through Roche-lobe overflow. In systems with a low 
white dwarf magnetic field, an accretion disk can form. The accretion disk is the origin for the 
strong emission lines characteristic of CV spectra not taken in outburst. In many cases, th e strong 



disk emission outshines the absorption features from the secondary star. IWarneii (|l995l ) gives a 
comprehensive review of CVs. 

Beyond observational difficulties, CVs (like most astronomical objects) present challenges in 
data interpretation. The determination of binary system parameters suffers from many potential 



^Based on Observations obtained at the MDM Observatory operated by Dartmouth College, Columbia University, 
The Ohio State University, and the University of Michigan 
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systematic effects. Typically, radial velocities from emission lines originating in the accretion disk 
are used to constrain the mass ratio of the stars, but these featu res do not always directly map 
the motion of the white dwarf (jWadd Il985l : iMarsh &: Hornd ll99Cl ) in part because of asymmetry 
in the brightness of the disk and line broadening due to Keplerian motion of the disk. When 
features from the secondary are observable, velocities may be calculated from its absorption lines. 
While the absorption lines generally have a larger velocity amplitude than the emission lines, 
inter pretation of absorption velocities is di fficult as the secondary star is by no means a "textbook" 



star ( Beuermann et al 



19981 : iKniggd 120061 ) . The donor star fills its Roche lobe and is therefore not 



spherical. Its surface is being irradiated on one side, as well. Despite these difficulties, eclipsing 
CVs are valuable because they provide the additional constraint on the inclination of the system. 

Our CV of interest, SPSS 154453+2553, was discovered i n the Sloan Digital Sky Survey 



(jAbazajian Sz Sloan Digital Sky Survey 



2008 



Szkody et al.ll2009l ). Spectra f rom SDSS show sin - 



gle peaked emission lines with a detectable contribution from an M-dwarf (ISzkodv et al.l |2009| ) 
consistent with our data. Photometry from SDSS gives a g magnitude of 16.60. 



2. Observations 

Our data consist of time-series photometry (see Section 2.1) and spectra (see Section 2.2). Our 
early spectra showed signs of a rotational disturbance in the emission lines and follow-up time-series 
photometry soon confirmed that the system eclipses. Radial velocity and eclipse observations taken 
in 2008 and 2009 were, when combined, sufficient to establish an unambiguous period of just over 
6 hours (see below). 



2.1. Photometry 

Time-series photometry was taken in June 2008, February 2009 and March 2009 at the MDM 
Observatory 1.3m telescope using a SITe 1024^ CCD and binned 2x2. Figure [1] shows differential 
V photometry taken during a June 2008 run at the MDM 1.3m telescope. Additional photometry, 
taken by Steve Brady with a 0.25 m robotic telescope in southern New Hampshire, was used along 
with the 1.3m data to determine an accurate ephemeris. 

The system showed small variations in brightness before and after eclipse with an eclipse depth 
of approximately one magnitude. There is no clear indication of the time of bright spot eclipse in 
our light curve. However, the feature at HJD 2454902.9 appears to be the white dwarf egress. The 
photometry and spectroscopy establish the ephemeris to be 

HJD mid-eclipse = 2, 454, 878.0062(15) + 0.251282(2)E (1) 

where E is the integer cycle count. The eclipse width, A(j) = 0.11 (~41 min.), was measured halfway 
down the steep sides of the light curve profile (Figured]), which corresponds to the brightest part 
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of the disk being covered and uncovered. 



2.2. Spectroscopy 



Spectra covering an entire orbit were taken in June 2008 with additional data taken in Au- 
gust 2008 and January 2009. These were obtained using the 2.4m Hiltner telescope at the MDM 
Observatory using the modular spectrograph and a SITe 2048^ CCD, which gave a resolution of 
3.5 A. Our coverage was from 4210-7520 A with vignetting toward the ends. Exposures were 
typically 300s in good seeing. We performed our wavelength calibration using comparison lamp 
spectra taken during twilight and shifts from the 5577 A night sky line. 

The average spectrum shows emission lines as well as weak but detectable absorption from the 
secondary M-dwarf (top line, Figure 2). The emission features appear single-peaked throughout 
the orbi t. To refine the spectral type of the secondary further, we subtracted a grid of M-dwarf 
spectra (lBoeshaaiill976l ) from our average spectrum and inspected the results by eye. At the best 
fit spectral type, we should be left with only continuum and emission features, while the M-dwarf 
contribution should ideally disappear. The best cancellation of the M-dwarf features occurs at type 
Ml. 5 lb 1.0 (see Figure [2]). By comparing the contribution from the M-dwarf in SDSS1544 to the 
template spectra used in the decomposition (with known ap parent rnagni tude), we can determine 



an apparent magnitude for the donor star of my = 19.0ib0.3. iKniggd (|2006l ) predicts a spectral type 
of K6.6 - earlier tha n our determi i iation . Our later spectral type determination implies that the 
secondary is evolved. iHoard et al.l (|200lh found a similar result in the system CM Phe. Although 
there was some ambiguity in their period determination, they d etermined a spectral ty pe of M2-5 
from a spectral deconvolution, while the predicted spectral type (ISmith &: Dhillonlll998l) is K 4-M1. 
The subject of late spectral types in CV donor stars is also discussed in iBaraffe fc Kdbl (l20od V For 
an orbital period close to 6 hours, their models predict possible spectral t ypes of K4-M2 depending 
primarily on q. In addition, two stars (XY Ari & LL Lyr) presented in Ixnigg^ ( 2006 )'s Table 2 
have orbital periods close to 6 hours with spectral type determinations of MO it 0.5 and M2.5 ± 1.5, 
respectively. 



Emission line y elocities from Ha were determined using a convolution method (ISchneider &: Young 



1980; IShafteiill983l ). In short, this method uses an antisymmetric function consisting of positive 
and negative Gaussians offset by an adjustable separation, a, which is convolved with the observed 
line profile. The zero of the convolution is taken as the line center. We adopted a va lue of a = 2080 



km s"^ for SDSS1544. To determine this value, we used the method developed bv IShafterj (jl983l ) 
to create a diagnostic diagram. This procedure allowed us to determine the maximum separation 
of the Gaussians for which the signal-to-noise ratio was sufficient to probe the line wings. 

Velocities from the absorption component of the sp ectrum were determined using the xcsao 
task in the rvsao package of IRAF feurtz &: Mink 1998 ). This was done by cross-correlating the 
spectrum between 5900 A and 6500 A with a composite M-dwarf spectrum with known radial 
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velocities (jMarcv et alJll987l ). 

Figure [3] shows the measured velocities from emission and absorption features along with their 
respective best-fit sine curves. These sine curves have the form v{t) = 7 + K sm[2'ir{t — To)/P] 
and Table 1 gives the fit parameters for each curve. Emission velocities affected by the rotational 
disturbance were excluded from the sinusoidal fit. If the emission lines were to represent the white 
dwarf motion, the eclipse should occur at (/> = 0.5 in the convention used here. However, there is a 
phase offset of Ac/) = 0.087ib0.007. This implies that the disk emission is somewhat asymmetric and 
these velocities are not faithfully representative of the motion of the white dwarf. The difference 
in the mean velocity between the emission and absorption curves also indicates that one is affected 
by something other than the stellar motion. Another interesting feature seen in Figure [3] is a clear 
rotational disturbance in the emission velocities seen as large deviations from the expected velocity. 
This is most obvious at eclipse ingress. A rotational disturbance occurs when the accretion disk 
is only partially eclipsed. As the secondary star moves in front of the rotating disk, we see only 
the redshifted side of the disk instead of the integrated light. During eclipse egress, the opposite 
phenomenon occurs and we observe only the blueshifted side. 



3. Analysis 

Normally, our observations would allow us to calculate the masses of the stars. However, this 
is complicated by the fact that the emission line velocities do not represent the true motion of the 
white dwarf. In addition, the absorption velocities show significant scatter with individual velocities 
having an uncertainty of 40 km s"^. This results in an error of 8 km s~^ in the absorption velocity 
amplitude and a greater than 10% in the mass calculation even if systematics were negligible. Even 
though the emission line velocities do not accurately represent the motion of the white dwarf, we 
will calculate the mass ratio for illustrative reasons. If we take the observed values as true indicators 
of the motion of the stars, we obtain q = K1/K2 = M2/M1 = 0.52. For a highly inclined object, 
this would imply a white dwarf mass of 0.45 Mq, which would give a companion mass of 0.23 Mq. 
This value is slightly lower than the typical value for an early M-dwarf star on the main sequence. 
We should be clear, though, that these values are likely not correct since the emission velocities 
do not track the white dwarf motion. However, this yields reasonable values for the masses of the 
stars and is a useful point of reference. Combining q = 0.52 and the eclipse half-width, we can 
constrain the inclination of the system to be z = 81° it 1°. 

We derive a distance to SDSS1544 through two different methods. The first combines our 
apparent magnitude from the spectral decomposition with an absolute ma gnitude found using the 



relationship between surface brightne s s and spectral type for normal stars (jBeuermann et al.lll999l ) 



(see Table [2]). From iBaraffe &: Kolbl (|2000l ) evolutionary models, we can find a R2 corresponding 
to our calculated mass for models with Po^b = 6.0 hr. Using their tabulated mass range for M2 of 
0.248-0.53 Mq, we find a radius range of 0.488-0.635 Rq. Adopting R2 = 0.55 ± O.OTPq, we can 
convert surface brightness to My and therefore find the corresponding My for our Ml. 5 dwarf. We 
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corrected for extinction using the infrared dust maps of lSchlegel et alj (|l998l ). These maps provide 
our best estimate for the amount of extinction due to dust. However, the method is imperfect as 
they give the amount of extinction to the edge of the Galaxy rather than to SDSS1544. Patchiness 
at scales below their 6-arcmin resolution would also affect the extinction estimate. Combining our 
values for my, My, and Ay , we c alcula te a distance of 800 it 180 pc. We may also apply the 
donor sequence derived by iKniggd (j2006l ) to calculate a distance (see, again, T a ble [2 ) • Using the 
K magnitude from the 2MASS survey and the absolute magnitude from iKniggel (j2006 )'s empirical 
donor sequence, we get a distance of 810 pc. 



4. Conclusion 



We find that SDSS1544 is an eclipsing CV with a period of 6.03 hr, and determine an ephemeris. 
The Ha emission velocities show a pronounced rotational disturbance. We also note that the 
determined emission velocities are not exactly out of phase with the absorption, which would point 
to asymmetric disk emission. The decomposed spectrum is fitted best with a Ml. 5 star, which 
indicates that the secondary is somewhat evolved. We also calculate the distance to SDSS1544 by 
two different methods. Both result in a distance of roughly 800 pc. 
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Table 1. Sinusoidal Fits 



Data Set 


Epoch (days) 


Period (days) 


K'^ (km s-i) 


7** (km s "'") 


(km s-i) 


N 


Ha emission 


54647.724(3) 


0.251282 


101(7) 


-44(5) 


34 


121 


Absorption 


54645.8182(15) 


0.251282 


193(8) 


11(6) 


40 


78 



'^K is the amplitude of the sinusoidal fit and 7 is the mean velocity. 
''Error in a single measure based on the scatter around the best fit 
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Eclipse Light Curve 

2009 March 12 




Time (HJD - 2454900.) 



Fig. 1.— The differential V light curve for SDSS 154453+2553. 



Table 2. Values used in Distance Calculations 



Method 


Absolute Magnitude 


Apparent Magnitude 


Extinction Distance (pc) 


Beuermann et al. (1999) 


9.29 ± 0.58(V) 


19.0 ± 0.3(V) 


0.2 ± 0.1 800 ± 180 


Knieee (2006) 


4.48(K) 


14.019 ± 0.065(K) 


— 810 
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I 

4500 5000 5500 6000 6500 7000 

Wavelength (A) 

Fig. 2. — A flux averaged spectrum in the rest frame of the secondary is plotted as the bottom 
spectrum. The average spectrum shows single peaked lines with a weak contribution from an M- 
dwarf secondary. The middle spectrum is the residual after subtracting a scaled Ml. 5 spectrum of 
Gliese 526 from our average spectrum, shown as the top spectrum. 
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Fig. 3. — Velocities from Ha are shown by plus signs while M-dwarf absorption is shown as filled 
squares. These arc shown with their best-fit sinusoids as the solid lines. The data are repeated for a 
second cycle for continuity. The velocities measured during the eclipse are not included in the sine 
fit. The velocities derived from the Ha emission are not exactly out of phase with the absorption 
velocities indicating asymmetric emission from the accretion disk. There is also a clear rotational 
disturbance seen in the Ha velocities. 



